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Characteristics of the high energy families (bundle of high energy e,7) and hadrons in the air-showers detected in 
the hybrid experiment together with emulsion chamber and AS-array at Mt.Chacaltaya are studied in detail by 
comparing with those of CORSIKA simulations using interaction models of QGSJET and EPOS. Because the 
atmospheric families and hadron component have more direct information of the nuclear interaction, correlations 
between atmospheric families and burst (hadron component of air-showers) accompanied to air-showers are more 
sensitive to the mechanism of the the cosmic-ray interactions. The burst size dependence of the family energy 
is compared with those of simulations. It is found that the family energy accompanied by the air-showers with 
the larger burst-size is systematically smaller than that expected in the simulated events. The experimental 
results can not be described simply by changing the chemical composition of primary cosmic-rays and this 
indicates that the x-distribution of secondary particles in cosmic-ray interactions becomes much steeper than 
that assumed in the simulation models. 



1. Introduction 



The main interests of the cosmic-ray study now shift 
to the highest energy region, Eq > 10 19 eV, using 
huge experimental apparatus to search for astrophys- 
ical sources and the acceleration or emission mecha- 
nism of those extreme high energy cosmic-rays. How- 
ever, half a century after its discovery, the "knee" in 
the cosmic-ray spectrum, a steepening of the energy 
spectrum at 10 15 ~ 10 16 eV, is still not well under- 
stood. Because of the low intensities, direct observa- 
tions of primary cosmic-rays in this energy range are 
still not possible and so various types of air-shower 
experiments at high mountains and also at ground 
level have been carried out in order to investigate the 
chemical composition of primary cosmic-rays in this 
energy region which gives important information on 
the physical origin of cosmic-rays [H-01- The exper- 
imental data in those indirect measurement are usu- 
ally interpreted by comparing with Monte Carlo sim- 
ulations assuming some models of cosmic-ray inter- 
actions. Many experimental groups claim that the 
fraction of heavy primaries increases rapidly beyond 
the "knee" region, e.g., the fraction of protons is esti- 
mated by the Tibet AS-7 group 8] to be as small as 
-10% of all particles for E = 10 15 - 10 16 eV. The 
results, however, depend on their assumed interaction 
model. For example, the EPOS model recently pro- 
posed [9|, Il0l | gives muon numbers much more than the 
QGSJET model. The events which can be interpreted 
due to heavy primary when we employ the QGSJET 
model as nuclear interactions are interpreted due to 
proton primaries when EPOS is used as the interac- 
tion model [ill ] . Thus the interpretations rely heavily 
on the Monte Carlo calculations. In fact, various ex- 
perimental groups give various data on chemical com- 
position in this energy region and the results are still 
very confusing. We should examine whether the over- 
all experimental data can be well interpreted by the 



assumed model before drawing a conclusion. 

Hybrid experiments operating simultaneously an 
air-shower array, a hadron calorimeter and an emul- 
sion chamber have been carried out at Mt. Cha- 
caltaya (5200m, Bolivia) (Ji], Yang-bajing (4300m, 
China) % H and Tien-Shan (3340m, Kazakhstan) Q] 
for studying cosmic-ray nuclear interaction in the en- 
ergy region around 10 15 — 10 17 eV. In the hybrid ex- 
periments, we can obtain the air-shower size, N e , from 
the air-shower array data, particle-density, rib, which 
are closely connected to the hadron component in the 
air-shower, from the hadron calorimeter (burst detec- 
tor) and energy and geometrical position of individ- 
ual high energy electromagnetic particles by the emul- 
sion chamber. Correlations between air-showers and 
accompanying families were studied so far in these 
three experiments by comparing experimental data 
and simulated data [IH, [l3j]. In the present paper 
we show some results obtained by studying correla- 
tions between the hadron component (data of hadron 
calorimeters) and families observed by emulsion cham- 
ber which are considered to be very sensitive to the 
mechanism of cosmic-ray interactions, using data of 
the Chacaltaya hybrid experiment. 



2. Short summary of the analysis on 
air-showers and accompanied families 

The shower-size, N e , dependence on the family en- 
ergy and on the lateral spread of the showers in the 
air-shower-triggered families were studied in the three 
hybrid experiments at high mountains. In Ref.[l3T] we 
have shown that the average family energy normal- 
ized by associating air-shower size, Ei? 7 /7V e , of the 
events with N e > 10 7 observed in these hybrid ex- 
periments agree more or less to those expected in the 
case of heavy-dominant composition of primary parti- 
cles, as shown in Figure [TJ though the difference in the 
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Figure 1: Depth dependence on the average family 
energy normalized by air-shower size for events with 
10 7 < N e < 10 8 . Solid lines are for proton-dominant 
composition and dotted lines for heavy-dominant 
compositionpji] . 
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Figure 2: Depth dependence on the average lateral 
spread of EAS-triggered families with 
l(f .V. 10 s 13. 



average value between the two chemical composition, 
proton-dominant and heavy-dominant, is not clear at 
the Tien-Shan altitude, especially in the case of the 
EPOS model. Some details of simulations are shown 
in Ref.[13(. The lateral spread of high energy show- 
ers in the families accompanied by the air-showers 
with N e > 10 7 was also studied 1 . The average lat- 
eral spread of showers in those families observed in the 
Chacaltaya experiment, shown in Figure [U is found to 
be smaller than that expected in the case of a heavy- 
dominant composition. The difference in the average 
lateral spread between the two chemical compositions 
is again not clear in the Tien-Shan data. The Cha- 



1 There is no official publication about the lateral spread of 
families accompanied by air-showers in Tibet AS-f experiments. 
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Figure 3: Configuration of 32 blocks of emulsion 
chambers and hadron calorimeters at the center of the 
Chacaltaya air-shower array. 
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Figure 4: The structure (side view) of one of the blocks 
of the emulsion chamber and the hadron calorimeter 
(burst detector). 



caltaya data show that the proton-dominant compo- 
sition is favorable to explain the small lateral spread 
of the families but the heavy-dominant composition is 
favorable to explain small family energy. Thus the in- 
crease of heavier composition of primary cosmic rays 
alone can not explain the general characteristics of 
air-shower-triggered families, contrary to the results 
of the Tibet group and others 0, H, Hi ■ 

The data of hadron calorimeters were also analyzed 
by the Chacaltaya and Tibet groups. The Tibet group 
concluded from the analysis that the experimental 
data were well explained by a heavy dominant com- 
position of primary particles but the Chacaltaya 
group concluded that the number of hadrons in the 
air-showers was less than expected [|| . 



3. Hybrid experiment at Mt.Chacaltaya 

The air-shower array covers a circular area within 
a radius about 50 m by 35 plastic scintillation detec- 
tors to measure the lateral distribution of the electron 
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Figure 5: Energy dependence of average number of 
charged particles arriving to the scintillator of the hadron 
calorimeter for 7r _ , proton, K~ , fi~ , e~ and 7 incidence. 




CORSIKA/QGSJET : 
proton E0=1.9e07 GeV, 
Ne=4.402e+06 
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Figure 6: Examples of burst data in the form of the 
lateral distribution of burst-density for experimental data 
(left) and for simulated data of proton-primary (right). 



4.2. High energy showers in emulsion 
chambers 



density of the air-showers. In the center of the air- 
shower array, 32 blocks of emulsion chambers (0.25 
m 2 each) are installed (see Figure |3|). Each block of 
the emulsion chamber consists of 30 lead plates each 
of 0.5 cm thick and 14 sensitive layers of X-ray film 
which are inserted after every 1 cm lead. The to- 
tal area of the emulsion chambers is 8 m 2 . Hadron 
calorimeters with plastic scintillator of 5cm thickness 
are installed underneath the respective blocks of the 
emulsion chamber (see Figure SJ . An iron support of 
2 cm thick is inserted between the emulsion chamber 
and the hadron calorimeter. Some detail of the Cha- 
caltaya hybrid experiment are described in Refs.P, 0] 



4. Simulations 



4.1. Air-showers and families 



For generating extensive air-showers and families 
we use the CORSIKA simulation code (version 6.735) 
[H employing the QGSJET model (QGSJETOlc) [Hf 
and the EPOS model (EPOS 1.60) 0, 03 for the 
cosmic-ray nuclear interaction. Primary particles of 
Eq > 10 15 eV are sampled respectively from the power 
law energy spectrum of integral power index —1.7, for 
pure protons and pure iron, and also from the energy 
spectrum of primary cosmic rays with proton dom- 
inant and a heavy dominant chemical composition. 
Some details of the chemical composition are shown 
in Table 1. The thinning energy is fixed to be 1 GeV. 
Shower size, N e , at the observation point is calcu- 
lated using the NKG option in the simulation. The 
air-shower center is randomly sampled within an area 
of ±2.5 m in the X and Y directions from the center 
of the hadron detectors (see Figured]). 



For high energy (e,7)-particles and hadrons of E > 
1 TcV in the atmospheric families arriving at each 
emulsion chamber, we calculate further the nuclear 
and electromagnetic cascade development inside the 
chamber taking into account the exact structure of 
the emulsion chamber. We use the QGSJET model 
for hadron-Pb interactions and a code formulated by 
Okamoto and Shibata for electromagnetic cascades 
(l9j . The electron number density under every 1cm 
Pb is transformed into spot darkness of the X-ray film. 
Then the energy of each shower is re-estimated from 
the shower transition on spot darkness by applying 
the procedure used in the experiments. 

4.3. Data of hadron calorimeters: 
Calculation of the burst-size 

Hadron calorimeters detect a bundle of charged par- 
ticles, which are produced in the emulsion chamber 
material by the hadron component in the air-shower 
through local nuclear interactions. Output from each 
unit of the hadron calorimeter is related to the en- 
ergy deposited in the scintillator, and is converted to 
a charged particle number using the average energy 
loss of a single muon in the scintillator. The number 
of charged particles per 50 cm x 50 cm, rib, is called 
the "burst density" hereafter. We use the GEANT4 
code [2(| for calculating the burst-density. We cal- 
culate the average number of charged particles 2 pro- 
duced in the emulsion chamber of 15 cmPb and ar- 
riving at the scintillator of the hadron calorimeter for 
hadrons (pions, proton, kaons), muons and high en- 
ergy e,7 in the air-shower, with 4 different energies 



2 Here we take into account the scintillator response of 
charged particles. Gamma-rays gives some energy deposit in 
the scintillator. Then the scintillator response of gamma-rays 
are also taken into account |2lH. 
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Table I Chemical composition of primary cosmic-rays and air-showers 

sampled primary particles 
proton dominant heavy dominant 



Eo (eV) protons He CNO heavy Fe 


protons He CNO heavy Fe 


10 15 - 10 16 42 % 16 % 16 % 14 % 12 % 
10 ie - 10 17 42 % 12 % 13 % 15 % 18 % 


17 % 10 % 18 % 15 % 40 % 
14 % 8 % 17 % 14 % 47 % 


air-showers accompanied by burst (C 
protons He CNO heavy Fe 


ORSIKA/QGS JET) 
protons He CNO heavy Fe 


10 6 < Ne < 10 7 57 % 18 % 11 % 9 % 5 % 

10 7 < Ne < 10 s 46 % 11 % 13 % 11 % 18 % 


31 % 13 % 11 % 13 % 26 % 
16 % 12 % 21 % 9 % 42 % 


air-showers accompanied by families of E_B 7 > 
protons He CNO heavy Fe 


10 TeV (CORSIKA/QGSJET) 
protons He CNO heavy Fe 


10 6 < Ne < 10 7 70 % 19 % 6 % 4 % 2 % 

10 7 < Afe < 10 s 50 % 9 % 14 % 7 % 20 % 


48 % 17 % 17 % 8 % 10 % 
23 % 10 % 19 % 3 % 45 % 



of 10 GeV, 100 GeV, 1 TeV and 10 TeV and 5 differ- 
ent zenith angles of arrival direction. Figure [5] shows 
an example of the energy dependence of average num- 
ber of charged particles which responds to scintillator 
for six different incident particles. The dependences 
are approximated by numerical functions and extrap- 
olate to higher or lower energy range of the incident 
particles. Applying these functions to every particle 
incident upon the emulsion chamber, we get the burst- 
density in each block of 32 hadron calorimeters. We 
define n™ ax as the largest burst-density among the 32 
blocks of the hadron calorimeters and En^ as the sum 
of burst-density of 32 blocks. In the following we pick 
up events which satisfy the following criteria; 

• iV e > 10 6 , 

• n r b nax > 10 4 , 

• Ras-Bs < lm, 

where Ras-Bs is the distance between the burst 
center and the air-shower center. 

The burst center is determined by the algorithm de- 
scribed in Ref.[2]. In the Chacaltaya data, 1,034 
events satisfy the above criteria in ~ 40 m 2 year expo- 
sure of hadron calorimeters. Among them, 73 events 
are accompanied by high energy atmospheric families 
of ££ 7 > 10 TeV (E min = 2 TeV). Figure 1 show ex- 
amples of experimental and simulated burst data. We 
can see the lateral distribution of the burst-density is 
well described by the power law function [2] . 

In Table 1, we show the fraction of proton, He, 
CNO, heavy and Fe components in the air-showers 
accompanied by families and also those accompanied 
by burst. In the shower-size region of N e > 10 7 , cor- 
responding to Eo » 10 16 — 10 17 eV, the fraction of 
each component is similar to those assumed in the pri- 
mary particles because almost all air-showers in this 
air-shower size region accompany families and bursts. 



5. Correlation between air-showers and 
bursts 

Figure [7] shows a scatter diagram between air- 
shower size, N e , and maximum burst density, n™ ax , 
of the event for the experimental data and for the 
simulated data of proton- and Fe-primaries. In the 
events of iron-primaries, n™ ax is more or less propor- 
tional to N e though n™ ax is weakly correlated to N e 
for the events of proton primaries. It is very natural 
because Fe-air-nucleus interactions are assumed to be 
the superposition of a number of low energy nucleon- 
air-nucleus collisions and so the fluctuation becomes 
small. But for proton-air-nucleus interactions, the po- 
sition of interactions and/or released energy at the 
interaction fluctuate widely event by event. The dis- 
tribution in the experimental data looks close to that 
in proton-primaries. Figure [8] shows distributions of 
n™ 1 /Ne for four different chemical composition of 
primary particles, pure proton, pure iron, proton- 
dominant and heavy dominant. The shape of the 
distribution for pure-iron primaries is very different 
from that for the others. There are almost no events 

with n^ ax /Ne >~ 0.02 (\og(n^ ax /Ne) > 1.6) in 

the iron-induced air-showers. On the contrary, a con- 
siderable number of events are found in this region 
of the distribution for proton-induced air-showers. 
There is no systematic difference in the shape for 
the other three chemical compositions, pure protons, 
proton-dominant and heavy dominant 3 , and also for 
the two different interaction models, QGSJET and 
EPOS. The experimental data are well described by 



3 Nearly half of the air-showers accompanied by bursts are 
due to protons and He-nuclei, even when heavy-dominant chem- 
ical composition is assumed in primary particles, as seen in Ta- 
ble 1. This is a reason why the difference in the shape of the 
distribution among these three chemical compostions is small. 
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Figure 7: Scatter diagram between air-shower size, N e , and maximum burst density, n™ ax , of the event. Shaded area in 
the simulated data is biased because the sampled primary energy is larger than 10 15 eV. 



T — ' — ' — ' — ' — T 




-4 -3 -2 -1 

log(n b (max)/Ne) 



Figure 8: Distribution of n™ ax /Ne. Circles are 
experimental data and lines are simulated data, solid 
lines : proton-primaries, dotted lines : Fe-primaries. 
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Figure 9: Distribution of En^/iVe. Circles are 
experimental data and lines are simulated data, solid 
lines : proton-primaries, dotted lines : Fe-primaries. 



the model calculation for these three chemical compo- 
sitions of primary particles. Figure [5] shows the dis- 
tribution of TiUb/Ne where En;, is a sum of over 32 
blocks of hadron calorimeters. Again we can see the 
experimental data are close to those expected for pure 
protons or mixed chemical composition of proton- 
dominant and heavy-dominant, though the number of 
events with smaller Snj, are less than expectation. An 
almost similar analysis was done by the Tibet group 
and they concluded that their data are well described 
by a heavy dominant composition (see Ref.Q). 



6. Correlation between bursts and 
families 

Figure [TU] shows a correlation diagram between 
n™ ax , maximum of burst size among 32 blocks of the 
event and accompanying family energy. The exper- 
imental data are compared with those of simulated 
data of proton-primaries and iron-primaries. As seen 
in the figure, the family energy is almost proportional 
to n™ ax in the simulated data irrespective of the pri- 
mary particles though the family energy of the events 
coming from iron-primaries are smaller than that from 
proton-primaries. In the figure we can see the fam- 
ily energy in the experimental data is systematically 
smaller than that of simulated data in the events with 
larger burst-density, n™ ax > 10 5 . Figure fTTI shows an 
average family energy versus n™ ax . It is clear that 
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Figure 11: Average family energy versus n™ ax in the 
burst-triggered families. 



the experimental data can not be explained simply by 
changing the chemical composition of primary parti- 
cles. 

One may argue that the smaller family energy in the 
experimental data is due to the systematic underesti- 
mation of the shower energy in the emulsion chamber. 
Figure [12] shows a comparison of the integral spectra 
of family energy observed by the present hybrid exper- 
iment and that by the emulsion chamber experiment 
of the Brazil- Japan collaboration}^ . As seen in the 
figure, the two spectra smoothly connect with each 
other. Thus we can conclude that the energy estima- 
tion of high energy particles in the emulsion chamber 
of the present hybrid experiment is not much different 
from that of the emulsion chamber experiment of the 
Brazil- Japan collaboration. 

One may also argue about an overestimation of the 
burst-density, especially beyond the region of n& > 
10 5 . But this possibility is also ruled out, because 
the distribution on n™ ax /N e is well described by the 
simulations as seen in Figure [8] 



: o families accompanied by EAS with Ne>1 A 6 : 

• all families observed in hybrid experiment 
- families observed by Brasil-Japan collaboration ■ 
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Figure 12: Comparison of the integral family energy. The 
solid line is for the emulsion chamber experiment of the 
Brazil- Japan Collaboration and solid circles are for 
air-shower triggered families and open circles are for all 
families observed in the present hybrid experiment. 



7. Discussion 

The previous analysis on the air-shower-triggered 
families shows that the average family energy of the 
experimental data is considerably smaller than that 
of simulations of proton-dominant primaries in the 
shower size region of Ne > 10 7 , as described in section 
2. In Figure[T3]we show a correlation diagram between 
family energy, ~EEj and associated air-shower size, N e , 
for the same events shown in Fig. 10. In the air-shower- 
size region of N e > 10 7 , the family energy in the ex- 
perimental data is systematically smaller than that 



expected for proton-primaries and those events look 
like coming from iron-primaries. 

The characteristics of the bursts accompanied by 
the proton-induced air-showers are very different from 
those accompanied by iron-induced ones, as shown 
in Figure [7j The considerable number of air-showers 
induced by proton-primaries accompany large burst- 
density which are not seen in the iron-induced air- 
showers, and the experimental data are close to those 
expected in case of proton-primaries. The contradic- 
tion of the above two arguments is well seen in the 
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Figure 13: Correlation diagram between air-shower size, 
N e , and accompanied lamily energy, S_B 7 , of the events 
shown in Figure [lOl 



correlation digram between bursts and families, shown 
in Figure [TU] and no model can describe the observed 
correlation. 

The spectra of high energy particles detected by the 
emulsion chambers are more sensitive to the mecha- 
nism of the particle production of the most forward 
region of the rapidity. The hadron calorimeter sup- 
plies the data of hadron component in the air-shower. 
The hadron component in the air-shower bears more 
direct information on the nuclear interaction than any 
other component such as electron and muon compo- 
nent in the air-shower, and the burst-size of hadron 
calorimeter gives a measure of interaction energy. 

Suppose the ^-distribution of the produced particles 
becomes steeper, the number of high energy particles 
detected in the emulsion chamber becomes small (i.e., 
detected family energy becomes smaller), because of 
the high threshold energy of the emulsion chambers. 
The hadron component detected by hadron calorime- 
ters, however, does not change much because of the 
lower detection threshold energy. Then the ration of 
Ei? 7 /rt™ a2: becomes smaller. The observed discrep- 
ancy between experimental data and simulated data 
can be explain in this way, that is, the experimental 
data indicate the ir-distribution of produced particles 
in the cosmic-ray nuclear interactions of Eq >~ 10 16 
eV is much steeper than that assumed in the models. 
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